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Abstract 

We consider a string model at one-loop related to a (T-model whose antisymmetric tensor field is con¬ 
structed as complex structure on the background manifold, specially on a manifold R x N where A is a 
complex manifold. As an example, we consider a homogeneous anisotropic (1 -I- 4)-dimensional cr-model 
where space part of the background is a 4-dimensional complex manifold. By solving the related one-loop 
/3-functions we obtain a static solution so that by reduction of this solution to (1-|-3)-dimension we obtain a 
static solution of Einstein equation where the matter sector is effectively interpreted as an inhomogeneous, 
anisotropic and barotropic matter satisfying all the energy conditions. Finally, the T-dual background of the 
solution is investigated and it is shown that the duality transformation and reduction processes commute 
with each other. 
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1 Introduction 

String theory is the most consistent theory for unifying the fundamental interactions, incorporating gravity. 
The cosmological implications of string theory (string cosmology) have recently received considerable attention 
mm- In this context, the string theory has been used in the study of the physical situation at very early 
universe [5]. 

The string theory in curved background, generally known as a conformal u-model, is concerned with the 
massless modes of the string including metric, a rank two antisymmetric tensor field and a dilaton field (j). 
In the cr-model approach, the conformal invariance condition is the vanishing of the /3-functions in all loops of 
the fields which may be interpreted as field equations of effective space-time action [HISJIS]. 

An interesting class of manifolds in the study of (super)string theory includes complex manifolds in which 
the set of solutions of string equations is based on complex non-Kahler manifolds [TJ |H1 [HI UHl HI]- In [H], 
Witten has constructed a topological a-model where the complex structure J of target space appeared in the 
Lagrangian in an interaction term like Kalb-Ramond field, through a supersymmetry. From another but related 
point of view, the flux H could be related to the almost complex structure by imposing supersymmetry at0 

dJ ^ *J7 [T?1 im [T5] . 

In this paper, we are interested in the idea of our recent paper to establish a physical interpretation for 
the complex structures m- We know that in general relativity the gravitational effects are characterized by 
a geometrical structure, so called metric, on a Riemannian manifold. But, on a complex manifold there is an 
additional geometrical structure, namely the complex structure J, which may be capable of carrying physical 
significance. In this regard, we have already presented a matter interpretation for the almost complex structure 
m- Here, following the same idea we will investigate a a-model where the role of antisymmetric H-field tensor 
is played by the complex structure (or more clearly by the fundamental antisymmetric Kahler form which is 
related to the complex structure with a hermitian metric). In this way, similar to the metric, the antisymmetric 
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tensor field in the context of cr-model will be a geometrical structure of the manifold. We will assume that the 
background space-time is a (1 -|- 4)-dimensional manifold such that the 4-dimensional space part is a complex 
manifold carrying a complex structure. Then, the B-field defined on the space-time is considered to be related 
with the complex structure of 4-dimensional part. In this way, the kinetic term of complex structure will 
appear in the effective action. At this level, we can obtain the solution of the lowest order of /3-functions in the 
non-critical dimension with a non-zero central charge deficit. 

The plane of the paper is as follows. In section 2, we start by collecting some preliminaries about string 
effective gravitational equations and Hermitian structure. In section 3, we introduce the special backgrounds 
of (1 -f 4)-dimensional cr-model consisting of a hermitian metric and a complex structure as i3-field. Then, the 
solution of /3-functions are investigated. In section 4, dimension reduction of the static (I -I- 4)-dimensional 
solutions to (1 -|-3)-dimensional space-time is performed. Then, we discuss the nature of derived solutions, with 
the particular anisotropic and inhomogeneous form of energy momentum tensor given by dilaton and B-field. 
In section 5, T-dual solutions are investigated. The paper ends with a conclusion. 


2 Review on cr-model equations and Hermitian structure 

Consider the two dimensional u-model in a curved background M, where its action is expressed by 

/ = ^ ( 1 ) 

such that the h denotes the metric of worldsheet E, gfio is the metric of U-dimensional space-time M, </) and 
Bjio are dilaton and antisymmetric tensor known as Kalb-Ramond field on M, respectively. The /t, v are D- 
dimensional target space-time indices and a,/3 are worldsheet indices and is the scalar curvature of the 
world sheet. 

In the leading-order of coupling a', the conditions of conformal (Weyl) invariance of cr-model ([T]), i.e. the 
one-loop /3-functions, are given by mum 

= ho - \ho - VaV4, (2) 

= y^ie^hop), ( 3 ) 

pW =-R+hH2^2VpV^$+{dp$f + A. ( 4 ) 

These equations are identical with the equations obtained by variation of the following effective action with 
respect to the fields g/ip, Bfj_c, and (j), respectively [HIUIS] 

S{gi,o, Bj^o, ^) = \J d^x^g e^R - - A). (5) 

Here, = UpppHo^o, where the field strength H of H-field is defined as follows 

HfiOp = dpBpp + dpBpp + dpBpp. ( 6 ) 

This action is so called string effective action at one-loop. 

In the non-critical theory, the constant A is related to the central charge deficit of the original theory, i.e. 
^ [T71 [T51 [H] . It is analogous to the non-vanishing cosmological constant in cosmology in standard 

theory of gravity uniiii]. In this case, the one-loop /3-functions requires some string scale curvature and this 
requires to consider higher-order corrections to the /3-functions ©-O- 

Vanishing of the one-loop /3-functions will provide us with a conformal cr-model at one-loop and it is at the 
ground of the so-called string cosmology at one-loop [3]. It is worth mentioning that this expresses the lowest 
order conformality and dose not indicate that the u-model is exactly-conformal. 
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The above formalism is in the string frame. Another useful frame, in which the string effective action ([S]) 
appears as the £)-dimensional Einstein-Hilbert action of general relativity, is the Einstein frame. The new 
metric is called Einstein metric and is obtained by a rescaling of the string metric as follows nmn] 

dfiv = e 9go- (7) 

where, the gfto denotes the Z?-dimensional Einstein metric whose Ricci scalar is R. So, in the the Einstein frame 
the form of the action will be 

S=IJ d^x^g {k - A). (8) 

In this frame, the equations = 0 in @ - O are transformed to the Einstein equations as follows 

= (9) 

where the energy-momentum tensors are given by 

- Ae"'^g^e^ , (10) 


2^(A) ^ 1 


- 

flKA 


:H^ggu 


( 11 ) 


2.1 cr-model equations at one-loop with Hermitian structure 

We consider a cr-model on a manifold M, at least one part of which is a complex manifold, so that the manifold 
carries another geometrical structure beside the metric, namely the complex structure J. Following the idea of 
giving a physical role to the complex structure, it is intriguing to investigate coupling of the complex structure 
to the string, included in the R-field interaction in the second term of the cr-model (HD- 

Before attempting to investigate such a model, let us collect some preliminaries about the complex structure. 
Let N be an even dimensional manifold, an almost complex Hermitian structure {g, J) on N with coordinates 
{a:*} consisting of a Reinmanian metric g = gijdx'^ 0 dx^ and a structure J = dx'‘ 0 dj satisfying (22] 

V J/ = -5l (12) 


gijdm'dn — ffr) 


(13) 


which are definitions of the almost complex structure and Hermitian metric, respectively. These two conditions 
introduce a Hermitian structure. Fundamental two form or Kahler two form H is defined by O = g[., J.) = 
^riijdx"^ A dxR whose components are given by = gikJj^. Then from dUD and (USD, we have 

^ij — ^ji. ( 74 ) 

Nijenhuis tensor of an almost Hermitian manifold is defined by [21] 

N/k = -Jr^djJ^r + Jr^dkJr - + Jr^rnJC- ( 15 ) 


An almost complex structure J on a manifold M is integrahle if and only if = 0; in this case the almost 
complex structure is called complex structure [22] . 

As mentioned above, we want to describe a cr-model on a special manifold where the H-field is related to 
the complex structure in a special case. For this propose, an appropriate part of the D-dimensional manifold 
must be a complex manifold equipped with a complex structure J. On the other hand, noting the fact that 
the components of J should be real themselves, through the equations (IT2]) and H3D, a diagonal metric with 
Minkowski signature is not capable of being a Hermitian metric associated with a complex structure. This 
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persuades us to consider a model on a manifold of the form M = R x N, where N is an even dimensional 
manifold which is demanded to be a complex manifold. The manifold N will describes the spatial part and 
therefore the signature problem will be removed. Summarizing, if the manifold M is denoted by hatted indices 
and N by the i, j,... indices, we are going to consider an antisymmetric B-field of the form 

B = ^Bjxodx^ A dx’' = ^ Vlijdx'' A dx\ Boi = 0, (16) 

where 7 is a coupling constant and fly = gjkJi^ is the Hermitian Kahler two form associated to gij. 

The basic ingredient of this model is therefore the dilaton and a complex structure in the form of B-field. 
This ansatz for i3-field has been already used in supersymmetry as a result of supersymmetry condition which 
imposes an ansatz on the dilaton field, too |12j . But, here we have used a different idea in choosing this ansatz. 
In fact, since we are interested in a strictly bosonic theory where the B-field is a geometric structure, for 
instance as a Kahler two form, hence the form of dilaton field is introduced merely by solving the /3-functions. 
From Einstein-frame point of view, the solution of one-loop /3-functions is an effective gravitational theory 
solution. However, by starting from a cr-model, not directly from the effective action, the origin of arising 
complex structure in the effective action becomes more clear. Furthermore, in the tr-model approach there 
is the advantageous of the T-duality transformation for constructing new classes of equivalent solution. As 
mentioned before, in fact the solution of lowest order of /3-functions dose not provide an exactly conformal 
(T-mo del. 


3 A static solution of the effective gravitational equations at one- 
loop in (1 + 4)-dimension 

This section is devoted to find a solution of the /3-function equations in a special case in which the 

H-field includes a complex structure. Before attempting to do this, let us remark some requirements. We are 
interested in an integrable complex structure, i.e. vanishing Nijenhus tensor, with VJ 7 ^ 0. 

In the particular dimension of (1-1-4), we take a general metric ansatz (in string frame) of the form 

ds^ = gjxodx^dx’' = —N {t) dt^ + gijdx^dx^ 

= -N (f) df + -^-^dr^ + (t) r^{de^ + sin^ 6dif^) + c^{t) dx^, (17) 

1 — kr^ 

where 

{/1} = {0,1,2,3,4}, {/} = {1,2,3,4}, 
and x^ = X being an extra space-like dimension. 

The 4-dimensional space will be required to be a complex manifold carrying a complex structure J, and in 
this context the gij is to be the associated Hermitian metric. So, the first step toward finding a solution for the 
/3-function equations (H])-® with the proposed H-field in (fTBl) is to determine the fiij (or equivalently the J^-^). 
Thus, we must solve the equations of complex structure (HU) and Hermitian metric condition (US along with 
the integrability condition, = 0 , in m- Solving them with the above metric ansatz gives the following 
Hermitian structure and H-field 

J = J/dx''®dj 

_ c( ) \/(l ^ dx® dr -\ - =dr 0 dx - ^—dO (Si — sin 9 d(fi (Si dg, (18) 

a[t) kr'^) sine 

B = {—kr^ + 1)”^ a (t) c{t) dr A dx + b (t)'^ sin9d9 A dip). (19) 

The predicted time dependent functions a(t), b{t) and c{t) in the above structure will be fixed by solving the 
/3-function equations. According to (ITOl) . the non-zero components of field strength tensor Hpgp can be found 
as follows 

Htrx = -7 \/(l - kr‘^)~^ {a{t)b {t) +a{t)b (t)J , (20) 
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HrSip = -2 7 (i) r sin 0, 


( 21 ) 


HtBip =-2 jc{t)c{t) sin 6. (22) 

The Weyl anomaly coefficients, the /3-function equations of the above metric and B-field with a consistent 
dilaton field (j) which is a function of {t, r) are given in the Appendix. 

After solving the /3-function equations, we get two set of solution where all the scale factor functions which 
are anticipated in space part of the metric m are fixed as constants with the following relations 

\/3 

a(t) = —, b{t) = V3qi, c(t) = gi, (23) 

where the qi is an arbitrary constant and the other function and parameters are given as following two type 
i: 

^ = 0, 7 = ^, ^(i) = —^ = Mr) + F{t) + (l>o, (24) 

where F{t) is an arbitrary function of time, A is a positive parameter and <^0 is constant part of dilaton and 
ii: 


^ = 0; <^ = ln(r)-f (/io, A = 0. (25) 

where the lapse function N (t) remains arbitrary in the second solution . 

The dilatons in these solutions are logarithmic function of r. This type of dilaton field specially appears 
in some black hole solution in string theory for example [5^. As mentioned before, if the A is a result of 
charge deficit in the theory, the higher order corrections of /3-functions must be considered. But, the first set 
of solutions which have a non zero A is not consistent with two-loop /3-functions. On the other hand, the A 
in critical dimensionality can be a result of flux or curvature in the remaining dimensions. But, this regime 
requires the cosmological constant (minimum of the potential) to be zero or negative. However, the first solution 
is not allowed again and the only accepted solution will be the second one with a zero A. 

Consequently, the explicit form of (1 -I- 4)-dimensional backgrounds achieved by considering a complex 
structure J in the space part of H-field, has the following form 

dsltringframe = -^(t) df + sffi^ Odif^) + (26) 


\/2 

B = sinddd f\dq) + dr /\ dx). 


(27) 


4 Dimension reduction to (1 + 3)-dimensions 


In this section, we are going to reduce the dimension of the (1 -|- 4) theory (with coordinates (t, r, 0, Lp, x)) to the 
(1 -|- 3)-dimensional space-time (with coordinates {t,r,0,ip)). By using the standard technique of Scherk and 
Schwarz [26] one can parameterize the 5-bein as follows [27] 


= 


q 


(28) 


Here, it is assumed that the fields are independent of the x coordinate, i.e. there is a killing vector qp defined 
by m 

q^dp = dx, (29) 
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where x is the flat version of the x coordinate and q = a/I f the fj, and v indices run over {0,1, 2, 3} then 
(1 + 4)-dimensional fields will be decomposed as following 


9xx 

VxxQ 1 ^XfJ, ^^1 


dxg. 

— 9xxQ = B^i/ -\- 



= gg.iy + Vxxq^A^A^, (^ = 0+iln9. 

(30) 

where the {gg,v, B^,y, </), A^, q} are (1 + 3)-dimensional fields. So, the expressions of them in 

dimensional fields are written as follows 

terms of (1+4)- 

9tiy 

- 9xg,gxi/ 1 2 ^ 1 1 " 1 

— gg.L' ^ ; 0 — 0 .In <7x9; ) 

9xx 4 


Bfj,i, 

_ D 1 9x[flBi/^x 

— “r ^ 5 — J^x^i 

9xx 


A, 

= 1^, 9=\9xx\^/^. 

9xx 

(31) 


Applying these on the metric, S-field and dilaton in (1^ . (I?7l) and (1^ we get the background fields as 

dsltringframe = (0 ^{3dr^ + dO^ + sin^ ddip^), (32) 


,^ = In (r) + - In ( 39 ^) + , 


(33) 


9 9 

Bev = —qi r sin0, 


B,= 


2 ’ 


(34) 


= q = i-- (35) 

oqi 

In turn, the (I + 3)-dimensional string effective action will contain the same terms as in ([S]) in which the 
(1 + 4)-dimensional fields are replaced by (1 + 3)-dimensional ones (I5^ - (I551) . 

In general, by a time rescaling dt' = N{t)dt in the metric (l32l) . the arbitrary N{t) function could be 
absorbed and the metric recast in the following form 

ds'itr^ngframe = dO"^ + sin^ Odif^). (36) 


We will use this rescaled metric in the following. 

In order to study the above result in the Einstein frame we may use a conformal transformation to obtain 
the metric in Einstein frame. So, in (1 + 3)-dimensions 



Einstein frame ^ string framei 

(37) 

or equivalently 


(38) 

whose Ricci scalar is 

33/4g-0o 

(39) 




Obviously, in the limit r —> oo we have R —> 0, which accounts for the asymptotically flatness. The 
corresponding Kretschmann scalar is obtained as 


K = 


31 \/3e^'^» 
36 qi ’ 


(40) 
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which shows an essential singularity at r — > 0. This singularity will be reconsidered in the following. 

Now, we consider the Einstein equations (jH]) for this space-time in order to construct an effective matter 
source. The Einstein equations are as follows 

+ rW), (41) 

where and are given by (10) and (11). So, the corresponding energy-momentum tensor (as a perfect 
fluid) is given bjo 

= P{r) + P.(r) + Pi{r) + e\e\), (42) 

where Pr{f) and Pi{r) are radial and lateral pressures. In this way one can obtain 

p{r) = Pi{r)=?.P2{r) = j^. (43) 

This energy-momentum tensor represents an anisotropic inhomogeneous barotropic matter with a positive qi. 
We know that the matter obeys the barotropic equations of state as follows 


P^{r) = Wrp{r), Pi{r) = wip{r), (44) 

such that the two positive constants are given by Wr- = 1 and wi = ^. 

Note that there are divergences in the pressures and Ricci scalar at r = 0 which may be interpreted by the 
presence of an effective charge introduced by the R-field and dilaton at r = cH. In other words, in this example 
the geometry is so deformed that there is effectively a charge at r = 0. The Pltvp is zero and the only non-zero 
component of field strength is the HrOtp, hence the nature of this charge is magnetic and not electric^- 
The total magnetic charge associated with H is calculated in Einstein-frame as follows 

M = f e^H = f sin(6*) dr Ad9 A dp. (45) 


Therefore, the magnetic charge per volume unit is g = 3^/'*(7^^^e'^”and hence the qi parameter in the metric is 
related to a magnetic charge. 

It is worth mentioning the energy conditions. One can show that the perfect fluid (43) simply satisfies all 
the energy conditions, including the null energy condition (NEC) p + Pi > 0, the strong energy condition (SEC) 
p -I- T,iPi >0, p -I- Pi > 0, the dominant energy condition (DEC) p > 0, \ p \>\ Pi \ and the weak energy 

condition (WEC) p > 0, p -I- Pi > 0. 


5 T-dual solutions 

T-duality (target space duality) [33] can be generally used in the cr-model context to generate a new class of 
solutions and background. The presented example at the previous section is independent of the two x and p 
coordinates, which will be regarded as isometry coordinates. In this section, we investigate the T-duality in 
(1 -I- 4) and (1 -I- 3)-dimensions. 

^The eT in non-coordinate basis are satisfying the relation 

= riaf,, = diag{-l, 1,1,1). 

In the metric II32II we have eT as = l,ed = e„® = -T—s. 

^In fact, this divergence is appeared in both S-field and dilaton energy momentum tensor. The divergence which is caused by 
the dilaton is originally the result of the particular form of dilaton filed as a logarithmic function of r ll33ll in this model. 

“^The Kalb-Ramond charge density vector which is visualized by electric charge is given by This charge is 

zero becouse of Hti/p = 0. 
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In (1 + 4)-dimensions, Buscher’s T-duality transformations [33] with respect to the isometry direction of x 
have the following fornll: 


QXX - ^ 5 Qx^ - ^ 

Qxx Qxx 


A _ 9x^ 

~ •) ^x^ — T 1 


B —B 


— 9^v> 


Qxx 

idxfl^l^x 

9xx 

{^9x^9xy Bx^Bxu) 


9xx 

^ = ^ + ln|fe|- (46) 

where g/ip, Bf^p and (j) are metric, antisymmetric tensor and dilaton fields of the following dual cr-model 
1 


I = 


Vh<fz{h°‘^gfipda.X^di3X^ + e^>^Bf,pd^X^dpX^ + 


(47) 


For the rescaled of model (1331) . (1371) after using transformation (1331) one can find the following dual metric 

string frame Q^pdxdx^ 


= —dt^ + <7j^(^dr^ + r'^dO^ + sin^ ddip"^ + 3 dx^ + ^^^drdx), 


and dual antisymmetric and dilaton fields 


_ fry 

B = qi sin (0) dO A dtp, 
^ = In (r) + i In (3gj + ^o- 


(48) 


(49) 


(50) 


The other isometry coordinate is the p coordinate. Similar to the (1431) . the T-dual transformation with 
respect to p in (l+4)-dimensions gives the following T-dual metric, antisymmetric and dilaton fields respectively 
as follows 


ds 


^ , = n-dT^dT^ 

string frame UjJj 


= —dt‘^ + ^qi^dr"^ + ^qi"‘r"‘d9‘‘ + 


1 


qi sim 9 


dp 


V2 

2sin9 


d9dp + -—y dx^, (51) 


3qi 


X V2 

B = -^dx A dr, 


(52) 

(53) 


^ = 31n (r) -I- 2 In (sin0) -|- + 2 In (gj). 

In (1 -I- 3)-dimension (32)-(35), the only isometry coordinate is p. We may apply two procedures here: we 
can find the T-dual backgrounds in (1 -|- 4)-dimensions and then dimensionally reduce the x coordinate; or, we 
may reduce the x coordinate and then find the T-dual solution with respect to p. In the second procedure, 
after dimension reduction of x coordinate, Buscher’s T-duality transformation with respect to p coordinate are 
given as follows 

- _ 1 - _ B,pm p _ Qipp 

9(p(p — 5 9^rn — 5 -^(pm — ? 

9pp 9pp 9pp 


Brn.n. — 


9mn — 9mn 


{.9pmBn(p 9'pnBm(p) 

1 

9pp 

{^9pxn9p>n B(p’imB(pn') 

5 

9p<p 

4> = (j)P\n\g^^\. 


(54) 


^We are using the notion of /3-function equations of ©-id where the dilaton is minus twice of the dilaton of the Ref. and 
m- Therefor, the T-duality transformation of dilaton of m with the form of 0 = 0 — 1/2 In \gxx\ is matched with this coefficient. 













Here, m and n indices indicate the {t,r,9} coordinates. Then using (I5^ - (I551) . the T-dual answers with respect 
to if coordinate are given as follows 


ds 


{string frame} 9 


dx^dx'' 


= -dt^ + + ^qlr^dO^ + ^ ^ ^ 

2 qfr^ sin (a) 


9 V2 

dv + . dedip, 

2 sin 0 


(55) 


(/) = 31n (r) + 2 In (sin0) + + 1/4 In (3) + 5/2 In (gj), (56) 

such that the T-dual 5-field is zero. On the other hand, if we apply the dimension reduction according to (131|) 
on the T-dual solutions dSU-dSSl), we will obtain the same results as (|55ll . (I56|) and similar to it, the reduced 
5-field would be zero, i.e. the dimension reduction procedure and T-duality transformation are consistent with 
each other. 


6 Conclusion 

We have developed a cr-model on a manifold RxN whose space part is a complex manifold. Specially, following 
the idea of giving a physical implication to the complex structure, we have considered the complex structure of 
the complex manifold as the antisymmetric tensor field of the cr-model in the form of 5 = ^ gkjJi^dx^ A dxK 
Therefor, since the metric field in cr-model is a geometrical structure, the role of antisymmetric filed is played 
by a geometrical structure on the manifold. Two static solutions are found for the one-loop /3-functions in the 
(1 -|- 4)-dimensions, one with a non-zero and positive A and another with a zero A. Considering the nature of 
A as a charge deficit in non-critical dimensions, the second solution is preferred. Then, its dimension reduction 
to (1 -|- 3)-dimensional space-time is performed and resulted in a static solution of Einstein equation with the 
energy-momentum tensor corresponding to an anisotropic inhomogeneous barotropic matter which satisfies all 
the energy conditions. From gravitational point of view, this is interpreted as a static solution of an effective 
gravity theory with a dilaton and antisymmetric 5-field representing matter sources. There is a divergence 
in the energy Ricci scalar, density and pressures at r = 0 which is interpreted as a natural result of the 
presence of an effective magnetic charge at r = 0. Finally, T-dual solutions are constructed in both (1 -(- 4) and 
(1 -|- 3)-dimensions. We have shown that in (1 -|- 3)-dimensional space-time the T-duality transformation and 
the dimension reduction procedures are compatible. 
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Appendix 

Here we give the non zero components of the /3-functions jdpc, related to the metric given by p7l) and the B-field 
given by ((H]), as follows @ 

/3qq = —{[(/ (t, r) -(- Ina(t) -I- 21n6(t) -f lnc(t)]3V (t) — N{t)[{\na{t) Inc(t))^ -I- 4 ln 6 (t)^]}a(t)^c (t)^ 6 (t)^ 

-1-2 c{t)'^N(t)a(t)b(t)[(l) (t, r) a (t) b(t) 2 a(t)b(t) b(t)d(t)] 2 c{t)a{t)'^c(t)b{t)'^N(t), (57) 

/3fi = 2 c(t )^6 {t) N{t){r^[N{t) {kr^ — l) </)" (t, r) -\- krcj)' it, r) N{t) -\- 1/2 a{tY^'^ -\- a (t) a{t)^ {t, r) a (t) d{t)] 

-f 2 [a(t)r^d(t)(i ( 7 ^ -I- l) Inc(t) -I- \nb{t) - ilnA^(t)) -|- ilnc(<)^r^ 7 ^a(t)^ -|- N{t) {k ( 7 ^ -|- l) - 7 ^)]}, 

(58) 

®The overdot and prime stand for diffrentiation ^ and respectively. 
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= r4>" [t, r) b (t) a(t) — a {t) (j)' (t, r)b(t)r — 2a (t) b(t) +2b (t) a (t) ( 7 ^ + l) , 


(59) 


P 22 = N{t)c{t)[b{t)^r{(kr^ — l) N{t)(l)' (t, r) + ln5(i) ^ (j) (t,r) + ln 6 (<) ( 27 ^ (^) 

+a (tf c(t)N(t)b (t) b(t)r^l-^lnN(t) + Ina(t) + Inc(t)] + 2 [(/c( 7 ^ + l)r^ - 7 ^ - 

(60) 

/3f3 = l/2sin2(0)/3|2, (61) 

/344 = a(t)^b(t)N(t){j^c(ty + 2c(t)c(t) + 2 c(t)c(t)[(j> (t,r) + ( 7 ^ + l) Ina(t) + 21n6(t) - l/21nA^(t)]} 

+c(t)^b(t)N (t) (62) 

= 2 N(t)a(t)b(t){c(t)^a(t)ln.a(t)^ — a{t)[c{t){\na{t) + \nc{t))]^^ (t, r) — c{t)'^d{t) — c (<) a{t)c{t)} 

+c{t)N{t)a{t)b (t) {a(i)[lnA^(t) — 41n6(t) + 2\nc(t)]c(t) + c (t) d(t){\nN(t) — 41n6(t))}, 

(63) 

Poi = {{kr^ - r) (j}' {t, r) - 4 + 2) (d(t)c (t) + a{t)c{t ))), (64) 

P 23 = c{t)b{t)^N (t) {r^a(i)^ [— 2 (0 + Ina(t) + lnc(t)]ln 6 (<) + rN{t) {kr^ — l) (j)' {t, r) 

—h\b{tYr'^a{tY + ln 6 (t)r^a [t)^ 0 {t, r) + N{t)} + b{t)N{t)c{t)b{t)a{t)'^r'^, (65) 


(3^ = —Ac{t)N{t)h{t)a (t) r‘^[b{t)a{t)c(t) + c(<)(d (t) b{t) + 2 a{t)b{t))] 

—4 [(fcr^ — 1 ) N{t)(j)" {t, r) + a{t)‘^4> {t, r) + l/2A^(t) (fcr^ — l) 0' (0 r)^ + 1/20(0 r)^ a(i)^]r^ 6 (t)^c(i)^A^(i) 
—4iV(t)^c(t)^ {?'&(i)^ (3fcr^ — 2 ) 0' (t,r) + [— 1 / 2 r^Aa (<)^ + /c ( 7 ^ + 1/3) r^ — 7 ^ — l] 6 (i)^ + l/3a(t)^} 
—Sb{t)c{t)'^ a{t)'^ N {t)r'^ {—1 / Ab{t)[—2hia(t) +lnA^(t) — 41n6(t) — 2 lnc(t)]0 (0 r) 

+1/8 6 (t) lna(i)^ 7 ^ + [1/4 6 (t) ( 7 ^ + 2) Inc(t) — 1/4 b{t)\nN{t) + 6 (t)]lna(i) 

—1/26(t)[—l/41nc(i)^7^ + [l/21nA^(t) — 2ln 6 (t)]lnc(i) +ln 6 (t)[(— 7 ^ — l) ln 6 (t) +lniV(i)]]}. ( 66 ) 
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